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Micro scale boiling behavior in the vicinity of graphite micro-fiber tips on the copper- 
graphite composite boiling surfaces is investigated. It is discovered that a large number of 
micro bubbles are formed first at the micro scratches and cavities on the copper matrix in 
pool boiling. In virtue of the non-wetting property of graphite, once the growing micro 
bubbles touch the graphite tips, the micro bubbles are sucked by the tips and merged into 
larger micro bubbles sitting on the tips. The micro bubbles grow rapidly and coalesce to 
form macro bubbles, each of which sitting on several tips. The growth processes of the micro 
and macro bubbles are analyzed and formulated followed by an analysis of bubble 
departure on the composite surfaces. Based on these analyses, the enhancement mechanism 
of the pool boiling heat transfer on the composite surfaces is clearly revealed. Experimental 
results of pool boiling heat transfer both for water and Freon-113 on the composite surfaces 
convincingly demonstrate the enhancement effects of the unique structure of Cu-Gr 
composite surfaces on boiling heat transfer. 
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= micro bubble diameter 
= macro bubble diameter 
= graphite fiber diameter 
= pressure 
= pressure difference 

= second principal radius of curvature of the vapor-liquid interface 
= first principal radius of curvature of the vapor-liquid interface 
= contracting rate of bubble throat radius 

= distance between coalescence plane and fiber tip 
= time 

= volume of microlayer in a micro bubble unit 
= vapor velocity 

= volume fraction of the graphite fiber in the composite material 
= height of the fiber tip plateau 
= average thickness of microlayer 
= maximum thickness of microlayer 
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I. Introduction 

A S is well known, the nucleate boiling is a greatly effective heat transfer mode. Although over years researchers 
have experimentally and theoretically studied boiling process, there are no theories or literature yet that can 
exactly explain the underlying heat transfer mechanisms. However, an accepted physics is that bubble nucleation, 
growth and departure control the heat and mass transport rate. The basic understanding of boiling phenomena on 
plain surfaces brought about ideas modifying the process through changing properties of the working fluid or/and 
the boiling surface to increase the boiling heat transfer coefficient. Thome 1 classified the enhancement techniques as 
passive and active techniques and considered that the passive techniques have the most practical importance. The 
passive techniques have been considered to include treated, rough, extended surfaces and surface tension devices. 
Fifteen years ago, a new passive technique, so called metal-graphite composite surfaces, was presented and studied. 2 
Since then a series of investigations have been performed to determine the enhancement effects of the composite 
surfaces on boiling heat transfer and their mechanism. 3 ' 7 It was disclosed that the nucleate pool boiling heat transfer 
in Freon-113 (a highly wetting liquids) on copper-graphite (Cu-Gr) composite surfaces is increased up to four to six 
times over on a pure copper surface. The critical heat flux (CHF) on metal-graphite heating surfaces have also been 
investigated 8, 9 and a correlation equation for predicting CHF has been suggested. 9 Liang and Yang 10 tested pentane 
(another highly wetting liquid) nucleate pool boiling performance on both Cu-Gr and aluminum-graphite (Al-Gr) 
composite surfaces with a fiber volume fraction of 50% and concluded that the presence of graphite fibers not only 
enhances overall boiling performance in developed region but also reduces hysteresis effects in the start up region 
which is important for the practical application of immersion cooling in the electronic industry. Chao 11 et al 
conducted experiments on nucleate pool boiling both in a highly wetting liquid, Freon-113, and a moderately 
wetting fluid, water, on Cu-Gr composite surfaces with different fiber volume fractions. The experimental results 
have shown that significant enhancement in pool boiling heat transfer performance in both highly wetting liquids 
and moderated wetting liquids like water occurs on Cu-Gr composite surfaces. 

The enhancement mechanisms of boiling on the composite surfaces are believed to be related with: 1) the higher 
thermal conductivity of micro-graphite fibers, 2) a large number of embryo cavities of appropriate size, and 3) the 
poorly-wetted property of graphite fiber tips. 10 In order to clarify the physics of the enhancement, it is important to 
have an understanding of the processes of bubble formation, growth and departure. 

For common solid surfaces, numerous models had been developed to predict bubble detachment diameters in 
pool boiling and were reviewed by Zeng et al. 12 Many experiments have showed that the bubble detachment had 
been implemented by a necking process. 13 17 Mitrovic 15 ' 16 analyzed the movement of three-phase-line (TPL) and the 
formation of bubble neck during growth of the bubble. He considered evaporation at TPL tended to create a convex 
interface, leading to establish the Laplace-pressure, which pushed the liquid radially outward to form the bubble 
neck. Kandlikar and Steinke 18 suggested that an intense evaporation at the TPL exerted a vapor recoil force pushing 
the interface into the liquid and consequently, a convex vapor-liquid interface was formed and developed to become 
the bubble neck. Although the bubble necking and broken-off processes were described by Mictrovic 15 and Mori and 
Baines, 17 but the dynamic process has not been properly analyzed. Recently, Zhang et al. 19 have analyzed the detail 
dynamic process of the bubble growth and necking, and subsequently given a determination method of the bubble 
departure diameter. As for the case of boiling on Cu-Gr composite surfaces, the bubble formation, growth, and 
departure dynamic processes have not been well investigated yet. 

II. Formation and Growth of Micro Bubbles 

Different from common metal surfaces, the metal-graphite composite surfaces are of their unique geometric and 
physical characteristics. The polished composite surfaces are smooth in macro scale but rough in micro scale. Yang 5 
investigated the surface structure through scanning electronic microscope (SEM). Direct top-view photomicrographs 
of a Cu-Gr composite specimen with volume fraction of 50 % graphite-fiber and fiber diameter of 8 ~ 10 pm clearly 
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display that the graphite fibers are embedded in the copper matrix perfectly without any gap or void, as shown in 
Fig. la. The photomicrographs also show that the fiber tips are rather smooth while the copper matrix surface is 
filled with a large number of microsized low-lying trenches and intermingled narrow grooves. A slightly tilted top- 
view of the composite surface shows that the fiber tips bulge out of the copper matrix, surviving the polishing due to 
their pliability and tenacity, and form individual plateaus distributing evenly on the composite surfaces, as shown in 
Fig. lb. 

in contradistinction to the intuition, the preliminary experimental observations have found that the micro bubbles 
originate from the valleys between the fiber tip plateaus instead of the fiber tips 5 . It is obviously attributed to that 
plenty of potential nucleate sites are created by the microsized trenches and intermingled narrow grooves on the 
copper matrix surface. On the other hand, although the fiber tips may have higher temperature than the copper 
matrix surface, the relative smooth surface of the fiber tip plateaus hardly provide initial nucleate sites. 

Once the growing micro bubbles touch the fiber tips, the micro bubbles tend to migrate to the tips from the 
copper matrix and firmly sit on the tips for the poor wetting character of the graphite and a large surface tension at 
the TPL. As the micro bubbles sitting on the fiber tips grow, the newly emerging embryo bubbles on the copper 
matrix would be swallowed by the micro bubbles at a touch, as shown in Fig. 2a. Both the rapid evaporation due to 
the higher temperature of the graphite fiber tips and the swallow of the emerging embryo micro bubbles speed up the 
growth of the micro bubbles on the fiber tips. The further growing micro bubbles inevitably coalesce with each 
other, as shown in Fig. 2b. The coalesced micro bubbles rapidly grow up and come into being micro mushrooms 
with their own microlayers forming a liquid network connected with bulk liquid under the coalesced micro bubbles, 
as shown in Figs. 2c and 2d. 

As analyzed by Zhang et al, 3 the micro bubbles coalesce at their maximum cross section. By reference to Fig. 2c, 
the micro bubble diameter at coalescence, D, is related to the fiber diameter, d, and the volume fraction of the 
graphite fiber in the composite material, a , which is identical with the area fraction of the fiber tips in the composite 
surface, by the following equation: 

nd 2 !4 = aD 2 (1) 

The height of coalescence plane of the micro bubbles is regarded as the maximum thickness of the microlayer 
§ m . It can be seen from Fig. 2b that 8 m = A+ S, where A is the height of the fiber tip plateau, estimated from Fig. lb 
at about d/3 0, S is the distance between coalescence plane and fiber tip, equaling J(p 1 2) 1 +(d/2) 2 = d^j ; r /(4a) - 1 / 2 • 



Figure 1. Photomicrographs (3,000 times) of a Cu-Gr composite surface with volume fraction of 50% 
graphite-fiber (8~10 pm diameter): a) direct top view; b) slightly tilted top view. 
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Coalescence plane 



c) d) 

Figure 2. Initiation and growth of micro bubbles: a) initiation and migration of micro bubbles; b) grow of 
micro bubbles firmly sitting on the fiber tips; c) coalescence of micro bubbles d) growth of coalesced micro 
bubbles and formation of micro-mushrooms. 


Therefore, 


8. = -J— -1+A 


2 V 4a 

Based on simple geometry, the volume of the microlayer in a micro bubble unit can be calculated by 


V,„ = D 2 S,„ --nS 


- D 2 +-d 2 +S 2 

4 4 


nd 2 
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Idealizing the liquid volume as a column of constant cross section, the average thickness of microlayer, 8, can be 
derived by V m /(D 2 -7rcr/4). As a result, the following expression can be obtained: 


8 = 


\ — nl 6-a/3 hr 1 

2(1 -a) \4a~ + 30 


(4) 


Obviously, both the maximum and average 
thicknesses of the microlayer are determined by 
the volume fraction of graphite fiber in the 
composite a and the fiber diameter d. As 
mentioned above, d = 8 ~ 10pm, simple 
calculation results for the maximum and 
average thicknesses of the microlayer at d = 
9pm are listed in Table 1. 


Table 1. Microlayer thicknesses at different values of a 


a 


0.1 

0.25 

0.5 

0.75 

5 m 

(pm) 

12.08 

6.89 

3.70 

1.28 

5 

(pm) 

6.10 

3.75 

2.41 

1.19 
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In boiling heat transfer, microlayers play an essential role. Thermal energy provided by the composite surface is 
first conducted through the microlayers and then delivered into the micro mushrooms as latent heat by evaporation 
at the microlayer surface. Because the microlayers are connected directly with the bulk liquid through liquid 
network, the evaporated liquid of the microlayers can immediately be replenished. It is the ample feed of liquid into 
the microlayers that enables the composite surfaces to have lower superheat at a given heat fluxes of boiling, i. e. 
enhancement of boiling heat transfer. 


TIT. Formation of a Macro Bubble 

The continually growing micro mushrooms inevitably coalesce and merge with each other to form a larger micro 
vapor mass, which is a macro bubble embryo, as shown in Fig. 3. When the bubble embryo is formed from several 
micro mushrooms, a macrolayer is created under the macro bubble embryo, connecting with the microlayers and 
bulk liquid through liquid network. As analyzed by Mitrovic 15 ' 16 and Kandlikar ls for a bubble growth, at the 
beginning of the bubble formation, the TPL of the bubble would move outwards to forms a convex vapor-liquid 
interface at the bubble bottom, and subsequently a bubble neck is created. The same process also occurs in the 
growth of a macro bubble embryo on the composite surfaces. As the macro bubble embryo grows, the TPL of the 
macrolayer moves towards liquid under the combination of high evaporation rate and capillary-pressure difference. 
Zhang et al' q analyzed this kind of TPL moving process in detail. In the case of a macro bubble embryo formed on 
the Cu-Gr composite surfaces, the moving TPL would stop at next circles lined with the corresponding neighbor 
graphite fiber tips for the poor wettability of the graphite fiber, and consequently, convex vapor-liquid interfaces are 
created at the macro bubble embryo bottom both the inner and the outer TPLs, as shown in Fig. 3b. The macro 
bubble embryo rapidly grows up and finally forms a macro bubble with a neck being of a throat radius of r t and the 
bottom radius of r h whose value is estimated at 0.1 ~ 0.2 mm, as shown in Fig. 4. 

IV. Departure of a Macro Bubble 

The strong evaporation of liquid from the micro and macro layers is firmly supported by the ample liquid supply 
through the liquid network at the bubble bottom, and creates larger and larger negative pressure at the throat, i. e. 
Ap th = pi, h - p v th > 0, where pi,h and p Vtt n are the static pressures of liquid side and vapor side at the throat, 
respectively. 

It should be noted that the static pressure difference varies along the bubble interface, as shown in Fig. 5. At the 
top of a departing bubble the static pressure difference AP T can be expressed as 4a/©,;, where ©/ is the bubble 
departure diameter, while the static pressure difference at the throat AP th is estimated by 



Coalescing micro mushrooms 



Figure 3. Formation of a macro mushroom from coalesced micro mushrooms: a) coalescence of micro 
mushrooms; b) formation of a macro bubble embryo. 
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A growing macro bubble 



Figure 4. Formation of a macro 
bubble 



Figure 5. Growth, necking, and departure of a macro bubble. 


A P„, = ( P,Wf h - p t rl )/ 2 - (cr / r th - cr / P„, ) (5) 

where p v and p\ are the density of vapor and liquid, respectively, W th is the vapor velocity at throat, r th is the 

contracting rate of the throat radius, i. e. liquid velocity at throat, cr is the surface tension, R,i, is the second principal 
curvature radius of the vapor-liquid interface at throat. 

As analyzed by Zhang et al, 19 the bubble necking and departure have three possible modes : 

a. If change rate of the static pressure difference at bubble throat is positive (d(Ap th )/dt > 0), the bubble neck 
would be broken off rapidly with a very short necking process. 

b. If d(Ap t) J/dt = 0, the bubble neck continually and gently contracts before the neck is broken off. 

c. For the case of d(Ap,i)/dt < 0, the bubble neck contracts to a certain diameter and then maintains it until the 
bubble grows big enough to be broken off at the throat. 

In mode a, since d(Ap lh )/dt > 0, the static pressure difference AP th becomes larger and larger, and consequently, 
the throat is rapidly compressed to zero and then the neck is broken off, while in mode b, the static pressure 
difference AP,i, maintains a constant till the buoyancy of the bubble exceeds the surface tension force at the throat 
during the throat contracting process and so the neck is broken off at a finite throat diameter. However, in the mode 
c, d(Ap lh )/dt < 0, the AP, h finally decreases to zero and the throat reduces to and then maintains a certain value till 
that the growing buoyancy exceeds the surface tension at the throat to break the neck off. For mode a, the ZIP,/, 
experiences a sudden change, from a certain value to an uncertain value when the throat diameter reduced to zero. 
As a result, the major portion of the bubble above the neck departures with a singular point at the broken throat 
(point A in the Fig. 5). For both modes b and c, a singular circle line A-A is created at the broken neck of the 
departing bubble, at which a theoretical infinity of AP, h occurs. This is the reason of occurrence of a liquid jet in the 
departure bubble. Of course, the corresponding singular point B in mode a or circle singular line B-B in modes b 
and c is formed at the residual vapor portion sitting on the composite surface, as shown in Fig. 5, but the strong 
evaporation from the liquid network at the residual bubble bottom withstands the action of the sudden change of 
AP lh , and at the same time rapidly fosters and develops a new bubble. 
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V. Enhancement of Boiling Heat Transfer 

It is because the unique structure of Cu-Gr composite surfaces that the strong evaporation in macro bubbles is 
supported by the ample liquid supply through the liquid network connecting the macrolayers and microlayers to bulk 
liquid, and consequently, much higher heat flux can be reached in the boiling heat transfer. 

Based on the force balance at neck of a departing bubble, Chao et al" have theoretically and experimentally 
revealed that higher vapor velocity produced by a larger local heat flux of boiling on the composite surface enables 
the bubbles to depart at smaller departure diameters, and consequently, results in an increased bubble departure 
frequency. This enhancement mechanism can successfully explain the experiment results of the enhanced boiling 
heat transfer both for F-l 13 and water. 2, 6 ' 1 1 


VI. Conclusions 

In pool boiling on Cu-Gr composite surfaces, micro bubbles originate first from the valleys between the fiber tip 
plateaus instead of the fiber tips on the graphite tips, and then migrate to the tips. The growing micro bubbles sitting 
on the fiber tips coalesce and merge each other to produce macro bubble embryos with a liquid network consisting 
of micro- and macro-layers. As an important parameter, the microlayer thickness is analyzed and formulated. 

It is the unique structure of Cu-Gr composite surfaces that creates a liquid network connected to bulk liquid, 
from whom an ample supply of liquid firmly supports the strong evaporation at the surfaces of micro- and macro- 
layers and fosters and develops macro bubbles. During the formation process of a macro bubble, a convex-concave 
interface is created at the bubble bottom under the combination of high evaporation rate and capillary-pressure 
difference, and then developed to a neck with a narrower and narrower throat. It is the strong evaporation that 
produces larger and larger vapor velocity and creates a great negative pressure at the bubble throat which accelerates 
the bubble departure with smaller departure diameters and higher departure frequency, enhancing the boiling heat 
transfer. 
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